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Abstract. We searched integrable 2D homogeneous polynomial potential with a 
polynomial first integral by using the so-called direct method of searching for first 
integrals. We proved that there exist no polynomial first integrals which are genuinely 
cubic or quartic in the momenta if the degree of homogeneous polynomial potentials 
is greater than 4. 



1. Introduction 

A Hamiltonian system with n degrees of freedom is integrable if the system admits n 
independent first integrals in involution (Liouville integrability) . It is a fundamental 
and important problem to investigate the integrability of Hamiltonian systems. Since 
the Hamiltonian itself is a first integral, the case of one degree of freedom is trivial. The 
simplest but non-trivial problem arises in the case of two degrees of freedom. For this 
case, the existence of an additional first integral guarantees the integrability. Let us 
consider a Hamiltonian system with two degrees of freedom, 



At present, there is no ultimate algorithm (necessary and sufficient conditions for 
integrability) to determine whether the system (|1 . 1| ) is integrable or not for a given 
potential V(x, y). 

The solution of the system is said to possess the Painleve property if it has no 
movable singular point other than poles. The Painleve property of the solution has 
been believed to be closely related to the integrability of the system, which is known as 
the Painleve conjecture (see, for example, JL2], [L3|, |(J). Although any rigorous relation 
between the Painleve property and the integrability has not been established, some new 
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integrable systems were detected || [OJ by postulating that the solution possesses the 
Painleve property (the Painleve test). 

For some Hamiltonian systems of the form it is possible to prove the non- 

integrability, i.e. the non-existence of an additional first integral. In the early 1980s, 
Ziglin p5| , p6| presented a non-integrability theorem and proved the non-integrability 
of some well-known Hamiltonian systems. Yoshida gave a criterion for non- 
integrability of Hamiltonian systems (|1 . 1|) with a homogeneous potential by using 
Ziglin's theorem [2^] . Recently, the differential Galois theory has become an important 
tool in attacking the problem of integrability (see, for example, [U). Morales- Ruiz 
and Ramis obtained a stronger necessary condition for integrability from their own 
theorem |TJ, [lTj based on the differential Galois theory (see also This necessary 

condition also justified the so-called weak-Painleve property |20| as a necessary 
condition for integrability for the first time |24[ . 

On the other hand, we have difficulty in proving the integrability. Even though a 
given system passes the Painleve test or satisfies the necessary condition for integrability, 
it remains unknown whether or not the system is actually integrable. In order to prove 
the integrability of the system, we have to present a first integral which is independent of 
the Hamiltonian. However, there is no general method to obtain an explicit expression of 
the desired first integral. Let $ be a first integral of the system ( |1 . 1| ) . Then the Poisson 
bracket of $ and H vanishes, which gives the following partial differential equation 
(PDE). 

dx dp x dp x dx dy dp y dp y dy ^ x dx dp x dx ^ y dy dp y dy 



In the present paper, we search polynomial solutions of the PDE ( |1.2j ), i.e. we 
assume that the first integral is a polynomial in (x,y,p x ,p y ) and that the potential 
is a polynomial in (x,y). The advantage of considering polynomials is that the PDE 
( |L2|) becomes an identity for (x,y,p x ,p y ) and that the problem reduces to completely 
algebraic one. In addition, we assume that the potential V(x, y) is a homogeneous 
polynomial, motivated by the following three points: (i) Any polynomial potential 
V(x, y) can be written in the form of the sum of homogeneous parts as 

V(x, y) = V min (x, y) H h V max (x, y), (1.3) 

where V m i n and V max are the lowest degree part and the highest degree part, respectively; 
and it can be shown jnj that if the system ( |1 . 1|) with the potential (|1.3|) admits a 
polynomial first integral, then the systems only with the lowest degree part and the 
highest degree part, given by 

H= \(pl+p 2 y ) +V min {x,y), H = -(j>l+p 2 y ) + V max (x,y) (1.4) 

both admit polynomial first integrals. Namely, in order for the system with a non- 
homogeneous potential to be integrable, each of the systems only with the highest degree 
part of the potential and the lowest degree part of the potential must be integrable. 
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(ii) As we will see in section [3[ the homogeneity of potentials assumes the weighted- 
homogeneity of first integrals, which simplifies the form of first integrals and reduces 
the complexities of computations. Indeed, the computations for non- homogeneous 
potentials are too complicated to deal with, (iii) The homogeneity of potentials plays an 
essential role to obtain the non-integrability criterions in fH| [23| , |24"f . For these reasons, 
we treat Hamiltonian systems ( |1 . 1| ) with a homogeneous polynomial potential of degree 

K 

k 

V(x, y) = V k = J2 <XjX k - j y j = a x k + a x x k ~ x y + ■■■ + a k y k . (1.5) 

j=0 

Let us here mention the rotational degrees of freedom. The integrability is preserved 
under rotations of coordinates, i.e. a potential obtained from an integrable potential by 
a rotation of coordinates is again integrable. We should identify such two potentials. In 
the present paper, we assume that a.\ — from the beginning, which partially removes 
the rotational degrees of freedom. 

Hietarinta || performed the direct search for integrable systems of the form (|1 . 1| ) , 
where the potential is a homogeneous polynomial of degree 5 or less and the additional 
first integral is a polynomial of order 4 or less in the momenta. The purpose of the 
present paper is to extend the degree of homogeneous polynomial potentials to an 
arbitrary positive integer in order to make a complete list of integrable systems in 
the range studied. 

We have to admit that our search does not cover the whole of integrable systems. 
In fact, there are some integrable systems with a rational or transcendental first 



integral || |I0 |. More generally, as seen in the approaches based on Ziglin's analysis 
or the differential Galois theory, integrability requires meromorphic first integrals. 
However, it would be almost impossible to single out all integrable cases without any 
restrictions on the class of the first integral and of the potential. It may be said 
that restrictive assumptions are indispensable for carrying out a thorough search for 
integrable systems. In this sense, our setting, in which the whole computations are 
tractable, is the fundamental one for a thorough search. 

This paper is organized as follows. The new result (theorem will be shown after a 
brief summary of the known facts in section |]. In section [| the statement of theorem |2| is 
amplified, followed by details of the computations in section f|. Finally, in section [|, we 
present a list of all integrable 2D homogeneous polynomial potentials with a polynomial 
first integral of order at most 4 in the momenta. 

2. The known facts and the new result 



The classification of first integrals up to quadratic in the momenta is well-known [ T0[] . 
The existence of a first integral linear in the momenta is related with symmetry, the 
invariance of the system under a transformation: the conservation laws of momentum 
and angular momentum correspond to the invariances by rotations and translations, 
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FT . First integrals 



respectively. These are particular cases of Noether's theorem 
quadratic in the momenta are exhausted by the Bertrand-Darboux theorem: 

Theorem 1 ( ||14|| ) The following three conditions are equivalent. 

(i) A Hamiltonian system of the form (|1.1|) possesses an additional first integral quadratic 
in the momenta. 

(ii) The potential function satisfies a linear partial differential equation called Darboux 
equation. 

(Hi) The system is separable in elliptic, polar, parabolic, or Cartesian coordinates. 

If we apply the Bertrand-Darboux theorem to the potential (|1 . 5| ) , then we obtain 
the following homogeneous polynomial potentials with an additional polynomial first 
integral quadratic (or linear) in the momenta. 

• separable in polar coordinates 

k/2 



V k = r k = (x* + y*) k/2 = £ 



m=0 



\ k -2m 2m 

y i 



m 



k 



even 



$ = yp x - xp y 
separable in parabolic coordinates Ul9 | 



V k = - 



r + x 



k+l 



-If 



r — x 



k+l 



[k/2] 

E 2 " 

m=0 



2m 



k — m 
m 



./>■ —2my2m 



X 



1 



^ = PyiVPx ~ xp y ) + -y Vfc_i 



(2.1a) 
(2.1b) 

(2.2a) 
(2.2b) 



• separable in Cartesian coordinates 

Vfc = Ax k + By k (2.3a) 
<2> = p 2 x + 2Ax k (including <P = p x for V k = y k ) (2.3b) 

Note that the potential separable in elliptic coordinates drops out because it can not 
be a homogeneous polynomial. It is also possible to obtain the above list by using the 
direct search (see f4|, |10|). 

There seems to be no special property of the system connected to the existence of 
first integrals of higher orders in the momenta. We can make polynomial first integrals 
of apparently higher orders in the momenta from the above polynomial first integrals 
(CT) , (ggg) , and O). For example, 

1 



# = (ypx - xp y ) , $ 



PyiVPx - XPy) + ~y Vfe_l 



<P = (p 2 x + 2Ax k 



k\2 



;2.4) 



are polynomial first integrals which are apparently quartic in the momenta for the 
potentials fl2.1aj ), ( |2.2aj ), ( |2.3a| ), respectively. On the other hand, Hall 0] and 
Grammaticos et al B found independently the potential of degree 3, 



lb 



(2.5a) 
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with an additional first integral 



1 .3 i 3 2 2 3 2 4 



$ = p y -^y PxPy + -^v p y --^y - -^y ■ (2.56) 



Ramani et al [19| found the potential of degree 3, 

V 3 = x 3 + -xy 2 + ^y 3 (2.6a) 
with an additional first integral 



o V^3 i a 1^2 /3 o V^3 i r>\ ( o /7T ■ 2 1 

$ = + - I 2~f ) PxP v + ( 3x y - V3«xy + -y 



'pi 



1 3 a , 2 4 1 5 , Sv 7 ^ 



+ + ^Vy 4 - ^xy 5 + (2.66) 



and the potential of degree 4, 



V, = z 4 + + iy 4 (2.7a) 

with an additional first integral 

* = Pj + ^ 4 p£ - 2xy 3 p^ + (^V + ^ 4 J p 2 + + + ^y 8 . (2.76) 

All of the three polynomial first integrals ( |2.56| ), ( |2.66| ), ( |2.76| ) are genuinely quartic in 
the momenta. Here, we mean by 'genuinely' that these first integrals cannot be reduced 
to first integrals of lower orders in the momenta. See also || for the discoveries of these 
three integrable cases. 

Now the following question arises. Are there any other potentials that admit 
a polynomial first integral which is genuinely quartic (or cubic) in the momenta? 
Hietarinta || searched integrable 2D homogeneous polynomial potentials of degree 5 
or less with a polynomial first integral of order 4 or less in the momenta by means of the 
so-called direct method and concluded that there was no integrable case other than the 
known ones in the range studied. We extended the degree of homogeneous polynomial 
potentials to an arbitrary positive integer. Specifically, we investigated the existence of 
polynomial first integrals which are cubic or quartic in the momenta for Hamiltonian 
systems of the form 

H = \{pl+P 2 y ) + a Q x k + a 2 x k ~ 2 y 2 + ■■■ + a k y k , (2.8) 

with the degree of the potential k > 3 (the cases for k — 1, 2 are obviously integrable). 
Note that the term of x k ~ 1 y in the potential vanishes for the removal of rotational 
degrees of freedom as mentioned in section |I[ As a result, we obtained the following 
theorem. 



Theorem 2 If k > 5, then the Hamiltonian system (|2.8| ) admits no polynomial first 
integrals which are genuinely cubic or quartic in the momenta. 

Therefore, the answer to the above question is "No." 
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3. Amplification of the statement of theorem |2] 

In this section, we give details of the statement of theorem Without loss of generality, 
we can assume that a polynomial first integral $ has the following properties. 

Property 1. A polynomial first integral $ is either even or odd in the momenta (p x ,p y )- 
Property 2. A polynomial first integral $ is weighted-homogeneous, i.e. $ satisfies 

^a 2 ^- 2 )*,^-^^ 7 ^ = a M $(x,y,p x ,p y ), (3.1) 

where a is an arbitrary constant and M is a constant called a weight. The property 
1 is due to the time reflection symmetry of the system. The property 2 is due to the 
scale-invariance of the system, which arises from the homogeneity of potentials. See 
Appendix A| for more details. It is easy to check that all the polynomial first integrals 



in section ^| satisfy the properties 1 and 2. 

3.1. Polynomial first integrals cubic in the momenta 

From the property 1, we can put a polynomial first integral which is cubic in the 
momenta in the form 

* = Aq(x, y)pl + A x (x, y)p 2 x p y + A 2 (x, y)p x p\ + A 3 (x, y)p 3 y 

+ B (x,y)p x + B 1 (x,y)p y . (3.2) 

If we regard the PDE ( |1.2|) as an identity for the momenta (p x ,p y ), then we have the 
following three sets of PDEs. 

' A 0x = 
A 0y + A lx = 

Aiy + A 2x = (3.3) 
A 2y + A 3x = 

A 3 y =0 

B 0x = 3A V X + A 1 V y 

B 0y + B lx = 2A 1 V X + 2A 2 V y (3.4) 
Bly = A 2 V X + 3A 3 V y 

B V X + B X V V = (3.5) 

where the subscripts x and y denote partial derivatives. From the property 2, the 
polynomials A{(x,y) are homogeneous of the same degree and so are the polynomials 
Bi(x,y). The PDEs (|3.3| ) have four homogeneous polynomial solutions, which classify 
the leading part of the first integral (|3.2|) into the following four cases. 



Case 1. # = a 3 (yp x - xp y ) 3 + B p x + B x p y 

Case 2. # = (a 2 p x + b 2 p y ){yp x - xp y f + B p x + B x p y 

Case 3. <P = [a x p\ + hp x p y + c^){yp x - xp y ) + B p x + B x p y 

Case 4. $ = a p 3 x + b p 2 x p y + cop x p 2 y + d pl + B p x + B x p y 
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Let us next consider the PDEs ( |3.4| ). We obtain the PDE for V(x, y), 
A 2 V XXX + {3A 3 - 2 A , )\ + (3A) - 2A 2 )V xyy + A 1 V yyy 
+ 2{A 2x - A ly ){V xx - V yy ) + 2(3A 0y - A lx - A 2y + 3A 3x )V xy 
+ (3A 0yy - 2A lxy + A 2xx )V x + (A lyy - 2A 2xy + 3A 3xx )V y = (3.6) 



Case 2 



by using 

d 2 (3A V x + AiVy) - d x d y {2A x V x + 2A 2 V y ) + d 2 x (A 2 V x + 3A 3 V y ) 
= d y B 0x — d x d y (B 0y + B lx ) + d x Bi y 

= 0. (3.7) 
For the above each case, the PDE ( |3.6| ) becomes 

Case 1 

a 3 {x 2 yV xxx + (2xy 2 - x z )V xxy + (y z - 2x 2 y)V xyy - xy 2 V yyy 

+ 8xyV xx + 8(y 2 - x 2 )V xy - 8xyV yy + 12yV x - 12xV x } = 0, (3.8a) 

{x V xxx ~\~ 4:xyV xxy -\- (3y 2x )V xyy 2xyV yyy 
+ 8xV xx + 16yV xy - 8xV yy + 12V X } + 
b 2 {-2xyV xxx + {3x 2 - 2y 2 )V xxy + kxyV xyy + y 2 V yyy 

- 8yV xx + \QxV xy + 8yV yy + \2V y } = 0, {3.8b) 

ai(2xV xxy + 3yV xyy - xV yyy + 8V xy ) + 

b\ ( xV xxx 2yV xxy -\- 2xV xyy -\- yV yyy 4:V XX -\- 4:V yy ) -\- 

ci(yV xxx - 3xV xxy - 2yV xyy - 8V xy ) = 0, (3.8c) 

Case 4 

a (3V xyy ) + b {-2V xxy + V yyy ) + c (V xxx - 2V xyy ) + ^(3^) = 0. (3.8d) 
The potential must satisfy the two PDEs (p.6|) and (|3.5|). They are transformed into 



Case 3 



identities for (x,y) by the substitution of the potential ( |1.5[ ). Therefore we finally 
obtain two sets of simultaneous algebraic equations for the coefficients of the potential 
and of the first integral. We searched their solutions for k > 3 under the condition 
that the leading part of the first integral does not vanish, i.e. 03 7^ 0, (a 2 ,b 2 ) 7^ (0,0), 
(ai,bi,Ci) 7^ (0,0,0), or (ao, b , Co, do) 7^ (0,0,0,0). Table |] shows the results for the 
four cases. 

3.2. Polynomial first integrals quartic in the momenta 

From the property 1, we can put a polynomial first integral which is quartic in the 
momenta in the form 

<P = A (x, y)p x + Ai(x, y)p\p y + A 2 (x, y)p 2 x p 2 y + A 3 (x, y)p x p\ + A 4 (x, y)p y 

+ B (x,y)p 2 x + B 1 (x,y)p x p y + B 2 (x,y)p 2 + C (x,y). (3.9) 
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Table 1. Integrable potentials with a first integral cubic in the momenta. 



Case 


Potential 


First integral 


1 

2 
3 

4 


Vk=r k = {x 2 +y 2 ) k/2^ 

v k = o 

V k =r k = (x 2 +y 2 ) fc / 2 , 
V k = x k 


k = even $ = (yp x — xp y ) 3 

k — even <J> — (yp x — xp y )H 

$ = b oPy (p 2 x + 2x k ) + d p y 



If we regard the PDE ( |i.2p as an identity for the momenta (p x ,p y ), then we obtain the 
following three sets of PDEs. 



.4 



Ox 







A(jy + A\ x 

A ly + A 
A 2y + A 
A 3v + A 



2x 



3x 



A 



l 3y 

iy 



lx 










(3.10) 







B 0x = AA V x + A 1 V y 

B 0y + B lx = 3AxV x + 2A 2 V y 

B ly + B 2x = 2A 2 V X + 3A 3 V y 



B 



a 
a 



2y 



Ox 



0y 



A 3 V X + AAiVy 

2B V X + B\V y 
B 1 V X + 2B 2 V y 



(3.11) 



(3.12) 



From the property 2, the polynomials Ai(x, y) are homogeneous of the same degree and 
so are the polynomials Bi(x,y) and Co(x, y). The PDEs ( p.lOj ) have five homogeneous 
polynomial solutions, which classify the leading part of the first integral ( p.9| ) into the 
following five cases. 

Case 1.0 = a 4 (yp x - xp y ) 4 + B p 2 x + B x p x p y + B 2 p 2 y + C { 
Case 2. <P = (a 3 p x + b 3 p y )(yp x - xp y f + B p 2 x 
Case 3. = (a 2 p 2 x + b 2 p x p y + c 2 p 2 y )(yp x - xp y ) 



C 



^ase z. <±> = ya 3 p x -f o 3 p y ){yp x - xp y y- -f n p~ + B x p x p y + B 2 p 2 y 
Case 3. # = (a 2 p 2 + 6 2 p^ + c 2 p 2 y )(yp x - xp^) 2 + 5op 2 + B x p x p y + ^p 2 + C 
Case 4. <P = (a^ + + cip^p 2 , + d x p\){yp x - xp y ) + EoP*. + Su?^ + B 2 p 2 y + C 

Case 5. # = a p* + &oj£p» + coPxP 2 , + doP*pJ + e Py + ^op 2 + B x p x p y + 5 2 p 2 + C 

Let us next consider the PDEs (|3.11|) . We obtain the PDE for V(x, y), 

A 3 V XXXX — 2{A 2 — 2A^)V xxxy + 3\A\ — A 3 )V xxyy — 2(2Ao — A 2 )V xyy y — AiV yyyy 

- (2A 2y - 3A 3x )V xxx + (QA ly - AA 2x - 3A 3x + 12A ix )V xxy 

- (12A 0y - 3A lx - AA 2y + 6A 3x )V xyy - (3A ly - 2A 2x )V yyy 
+ (3A lyy - AA 2xy + 3A 3xx ){V xx - V % 



-(12 A, 



Oyy 



yyj 

6A lxy + 2A 2xx — 2A 2yy + QA 3xy — 12A ixx )V xy 



{AAo yyy 3A\ xyy -\- 2A 2xxy A 3xxx ^V x 
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~~ i^lyyy — 2A 2xyy + 3A 3xxy — ^A^ xxx )V y = (3.13) 

by using 

d*(4A V x + A 1 V y ) - d x d 2 (3A 1 V x + 2A 2 V y ) + d 2 x d y (2A 2 V x + 3A 3 V y ) - d 3 x (A 3 V x + AA 4 V y ) 
= d y B 0x — d x d y (B 0y + Bi x ) + d x dy(Bi y + B 2x ) — d x B 2y 
= 0. (3.14) 

For the above each case, the PDE (|3.13|) becomes 

Case 1 

a>^{x yV xxxx {x 3x y )V XXX y 3(x y xy )V xxyy + {y 3x y )V xyy y xy V yyyy 
+ 15x 2 yV xxx - 15(a; 3 - 2xy 2 )V xxy + 15(y 3 - 2x 2 y)V xyy - 15xy 2 V yyy 
+ 60xyV xx - 60(x 2 - y 2 )V xy - 60xyV yy + 60yV x - 60xV y } = 0, (3.15a) 

Case 2 

&3{% V xxxx 4" 6x yVxxxy 3(x 3xy )V xxy y 2(3x y 2y )V xyy y 3xy Vyyyy 
+ 15x 2 V xxx + 60xyV xxy - 15(2x 2 - 3y 2 )V xyy - 30xyV yyy 
+ 60x14^ + !20yV xy - mxVyy + 60V X } + 

b 3 {-3x 2 yV xxxx + 2(2x 3 - 3xy 2 )V xxxy + 3(3x 2 y - y 3 )V xxyy + 6xy 2 V xyyy + y 3 V yyyy 

- 30xyV xxx + 15 (3x 2 - 2y 2 )V xxy + 60xyV xyy + 15y 2 V yyy 

- 60yV xx - 120xV xy + 60yV yy + 60V y } = 0, (3.156) 

Case 3 

a 2 {2x 2 

VxXXy ^•^yV XX yy 2{X 2y )V X yyy 2XyVyyyy 

+ 20xV xxy + 30yV xyy - WxV yyy + A0V xy } + 
b%{ V xxxx Axy V xxxy ~\~ 3(x y )V xxy y ~\~ 4:XyV X yyy -\- y V y y yy 

- 10xV xxx - 20yV xxy + 20xV xyy + WyV yyy - 20V XX + 20V yy } + 
c 2 {2xyV xxxx - 2(2x 2 - y 2 )V xxxy - 6xyV xxyy - 2y 2 V x 



mm 



Case 4 



+ WyV xxx + 3QxV xxy - 20yV xyy - 40V xy } = 0, (3.15c) 

&\{^>%Vxxyy ^yV X yyy %Vyyyy H~~ l^V^y^) 

^1 ( ^^^xxxy ^IJ^xxyy ~t~ ^^^xyyy ~t~ V^yyyy lO^crry ^^yyy) 

Cl^XVjr.jr.jr.jr. ~\~ ^"{jV XXX y ^XV XX yy ^]jV X yyy ~\~ 5V XXX -LOV^yy) ~\~ 

d\{-yV xxxx + 4xV xxxy + 3yV xxyy + lW xxy ) = 0, (3.15c?) 



Case 5 



(^0^V X yyy) ~\~ ^o( 3V XXy y -\~ Vyyyy) -\~ C()(2V XXX y 2V X yyy) 

+ d (-V xxxx + 3V xxyy ) + e (-W xxxy ) = 0. (3.15e) 
We also obtain the PDEs for V(x,y), 

Bi(V xx - Vyy) + 2{B 2 - B )V xy + (B lx - 2B 0y )V x + {2B 2x - B ly )V y = (3.16) 
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Table 2. Integrable potentials with a first integral quartic in the momenta. 



Case Potential First integral 



1 


v k 




= (x 2 + 


y 2 ) k/2 , 


k = even 


<2> = 


{yPx ~ 2:p a ) 4 






2 


v k 


= 
















3 


v k 




= (x 2 + 


y 2 ) k/2 , 


k = even 


= 


(yp x - xp y ) 2 H 








v k 


1 

r 




-' + (-! 


k ^r— 


3> = 


(Py(yPx -xp y )-\ 




-0 2 


4 


v h 


_ 1 

r 


^r+x^ k ^ 


- X + (-l 


k ^ r—x j 


<2> = 


(Py(yPx ~XPy)-\ 


- b 2 v k - 




■5 


v k 


= Ax k + By k 






<Z> = 


a (pl + 2Ax k ) 2 - 


f e (^ - 


Y2By k ) 2 



from ( |3.12| ) by using 

9„(2S ^ + B X V V ) - d x (B 1 V x + 2S 2 V y ) = d y C 0x - d x C 0y = 0. (3.17) 

The potential must satisfy the two PDEs (|3.13|) and fl3.16|) . They are transformed into 
identities for (x,y) by the substitution of the potential (|1.5|) . Therefore, we finally 
obtain two sets of simultaneous algebraic equations for the coefficients of the potential 
and of the first integral. We searched their solutions for k > 3 under the condition 
that the leading part of the first integral does not vanish, i.e. a 4 ^ 0, (a 3 ,6 3 ) ^ (0,0), 
(a 2 ,& 2 ,c 2 ) ^ (0,0,0), (ai,6i,ci,di) ^ (0,0,0,0), or (a , b , c , d , e ) ^ (0,0,0,0,0). 
Table ^| shows the results obtained from the five cases, which hold for k > 3 except that 
there exist three exceptional potentials, (|2.5o|) , (|2.6o|) , (|2.7o|) , in the case 5 for k = 3,4. 



All the potentials in tables [I] and ^| are separable ones given in section |2] and the 
corresponding first integrals are apparently cubic and quartic in the momenta, i.e. there 
are no new integrable cases. This is what the statement of theorem |2| means. 

4. Details of the computations 



All of the three polynomial first integrals ( |2.5£|) , (|2.66|) , ( |2.7£|) , which are genuinely 
quartic in the momenta, fall into the case 5. So, it is quite natural to expect that new 
integrable potentials, if any, would have first integrals which belong to the case 5. For 
this reason, we take the case 5 as an example to show the details of the computations 
performed in this study. We first exclude the square of the Hamiltonian by considering 
<E> — 2cqH 2 instead of $ itself, i.e. we put Cq = from the beginning. The PDEs (|3.11| ) 
for this case become 



B Qx = 4a V x + b V y 
Bo y + B\ x = 3boV x 
B ly 
B, 



>2y 



B 2x = 3d V y 
d V x + 4:e V y 



(4.1) 
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We obtain the expressions of B , B%, B 2 by integrating ( |4.1| ) as follows (note that a\ = 0). 
fe— 1 

4(fe - j)aja Q + (j + l)a j+1 b k _ 



3=0 

Bi = rix k - kr 2 x k ~ 1 y 

k 

+E 



a;* V +r y*, 



fe - j 



(4.2) 



3=2 



(fc-j + l)(fc-j + 2) \ 4(fe-j + l) 

OQ!,- — ~ «j_2 f «o ; «j-ieo 



i0 - 1) 



x^V, (4.3) 



B 2 = r 2 x k + 

3=1 

where 



(fe ~ J + ljoj-jjo + 4jqje fc.j- j- 



r 



r 2 



3(fe - l)(fe - 2)a k ^b + {6a fe _ 3 - 3(fe - l)(fe - 2)a fc _i}d + 8(fe - 2)a fc _ 2 e 



fe(fe- l)(fe-2) 



8(fe — 2)a 2 a + 6(^360 
fe(fe-l)(fe-2) ' 
{3fe(fe — l)«o — 2a 2 }& 



(4.4) 

(4.5) 
(4.6) 
(4.7) 



fe(fe-l) 

The PDE ( gig ), or ( fUgg) becomes 

^OjQ^xyyy ^o{3V^xyy ^yyyy) do(V X xxx 3V X xyy^) A^O^xxxy 0. (4.8) 

From ( |4.8| ) and Q3.16j ) with B ,Bi,B 2 given above, we obtain two sets of simultaneous 
algebraic equations for the coefficients of the potential, aj, and of the leading part of 
the first integral, a , b , d , e . The simultaneous algebraic equations obtained from ( |4.8| ) 
consist of the following fe — 3 equations. 

±(j -l)j(j + l)(k-j -l)a j+1 a 

+ J(J ~ 1){0" + 1)0" + 2 )«3+2 - 3(fe ~j)(k - j - l)a 3 }b 

+ (fe - j)(fe - j - i){3j(j - 1)«3 - (fe - J + i)(fe - J + 2K_ 2 H 

- 4(j - l)(fe - j - l)(fe - j)(k - j + l)«,_ ie o = (j = 2, 3, . . . , fe - 2) (4.9) 

which can be regarded as a recurrence relation for aj. The simultaneous algebraic 
equations obtained from ( |3.16 ) consist of 2fe — 1 equations of the form 

M u bo = 
M 21 b + M 22 a = 
M 31 b + M 32 a + M 33 d = 
M 41 6 + M 42 a + M 43 d + M 44 e = 
M 51 6 + M 52 a + M 53 d + M 54 e = 
M 61 b + M 62 a + M 63 rf + M 64 e = 



(4.10) 



M2fc_i,i&o + M 2 fe_i )2 a + M 2fe _i )3 d + M 2A ._ M e = 
where the first four equations are given by the following. 
{3fe(2fe - 1)oq - 2ct 2 }{fe(fe - 1)oq - 2ct 2 } 



Mi 



fe(fe-l) 
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M 21 = x — + 

M 22 = 



k-2 k(k-l)(k-2) 
16(3fc - 2){k(k - l)a - 2a 2 }a 2 



k(k-l) 

2(17k 2 - 28k + 6)a 2 18(7 k - 6) a* 
M 31 = 6(k - l)(2k - 3)a a 2 - — rjz tt — — + 



k(k-l) k(k-l)(k-2) 
12(7k 2 - 22k + 18)a tt4 48(2Jfe 2 - 4A; + 3)a 2 a 4 



(fc - 2)(fc - 3) fe(As — 1)(A; — 2)(A; — 3) ' 

12(2A;-3)(3A;-4)aoa3 2A(2k - 3)(5Jfe - 4)a 2 a 3 
32 ~ A^2 + k(k- l)(k-2) ' 

M 33 = 4k(2k - 3)a a 2 , 

m ad- ovoi. i\ 4(29fc 2 - 44fc + 6)a 2 a 3 , 48(7*; - 12)a 3 a 4 
= 6(* - 2)(2Jb - 3)a a 3 ^-^ + _ 2)(fc _ 3) 

20(7fc 2 - 31Jfe + 36)a tt5 40(5£; 2 - Ilk + 12)a 2 a 5 
(fc-3)(fc-4) + k(k- l)(fc-3)(fc-4) ' 
16(fc - 2)(3Jfc -4)a 2 144a 2 ; 96(£; - 2) 2 a a 4 384(A; - 2)a 2 a 4 
42 ~ k~^l + ~k~^2 k^3 + k(k-3) ' 

M 43 = 12k(k - 2)a a 3 , M 44 = 32(A; - 2)a 2 

Now, what we have to do is to find solutions of ( f4.9|) and (|4.10|) under the condition 
(bo, a , d , e ) ^ (0, 0, 0, 0). Applying this condition to the first four equations of (|4.10|) , 
we see that the product MnM 22 M 33 M 44 must vanish. We therefore have the following 
possibilities for the relation between a and a 2 . 

(i) a = (ii) a 2 = (iii) a 2 = a (iv) a 2 = a (4.11) 

Once the relation between a and a 2 is given, the computations to follow are quite 
straightforward. Let us move on to the details of each case. 



(i) When ao = 0, the first and the second equations of fl4.10Q become 



4a 2 

J o 



2 b n = 



k(k — 1) , , 

12(7fc - 6)a 2 a 3 32(3*; - 2)a\ 1 ' 

k(k-l)(k-2) b ° + k(k-l) a ° ~ ° 

If we assume that a 2 ^ then we see that bo = ao = 0, and then we can show from (|4.9| ) 
with j = 2,3 that d = e = 0. This contradicts the condition (b , a , d , e ) ^ (0, 0, 0, 0). 
Therefore, a 2 must vanish. Then the third equation of ( |4.10|) becomes 

,^:„ 6 > Q V = 0. (4.13) 
fc(fc-l)(fc-2) V ; 

Let us assume that a 3 ^ 0. Then, b = and we can show from ( |4.9j ) with j = 2,3,4 that 
a = d = eo = 0. This contradicts the condition (bo, ao, do, eo) ^ (0, 0, 0, 0). Therefore, 
a 3 must vanish. Let us now suppose that we have proved that ao = «i = • • • = a? = 
(j > 3). Then it follows from ( |4.9| ) that atj+x must vanish, up to j = k — 4, under the 
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condition (b , a , d , e ) 7^ (0,0,0,0). Now, we have a = a\ = ■ ■ ■ = a k -3 = 0. If we 
assume that a k -2 7^ then we have bo = a o — do = from ( |4.9| ) with j = k— 4, k— 3, k—2. 
Then the (2k — 4)-th equation of ( 4.10 ) becomes 



16(fc-2)(3fc-4)a 



2 

fc-2 



, | e = (4.14) 

which shows that eo = 0. This contradicts the condition (60, ao, do, eo) 7^ (0,0,0,0). 
Therefore, 0^-2 must vanish. Let us assume that a k -i 7^ 0. Then we have 60 = «o = 
from ((L9|) with j = Jfe - 3, k - 2. The (2Jfe - 3)-rd and the (2k - 2)-nd equations of (|4TT0| ) 
become 

' -3(jfe-l)(2jfe-3)a^_ 1 d = 

fitt. 1V91- 1^ * , 8(fe-l)(3fc-lK_x n (4-15) 
-6(fc - l)(2fc - l)a fc _ia fc cto H e = 

which show that do = eo = 0. This contradicts the condition (bo, ao, do, eo) 7^ (0, 0, 0, 0). 
Therefore, a k -i must vanish. It has shown that ao = ct\ = ■ ■ ■ = otk-i — 0. Then the 
(2k — 3)-rd and the (2k — l)-st equations of Q4.10| ) become 

h 2 (k-l)(2k-Z)a 2 k bo = () (416) 
-3A;(2A; - l)a 2 k d = 

The only possible solution is bo = do = with a k 7^ 0. Therefore, we obtain 

V k = y k . (4.17) 

If we assume that ao = for the other three cases (ii), (iii), (iv), then we only obtain 
the potential ( J4.17[ ). Therefore, we assume that a 7^ for the cases (ii), (iii), (iv). Let 
us put a = 1. 

(ii) When a 2 = and ao = 1, we can see from the first equation of (|4.1C| ) that 60 = 0. 
Then the third equation of fl4.10|) becomes 

12(2k-3)(3k-A)a 3 
^— 2 a o = 0, (4.18) 

which shows that o^ao = 0. Then the recurrence relation ( |4.9| ) gives 

-Jfe(jfc-l)(A;-2)(fc-3)do = (j = 2) (4.19) 

from which we obtain do = 0. Let us now suppose that we have proved that 
olx — «2 = • • ■ = Oij — (j > 2). Then it follows from (|4.9| ) that aij+i must vanish, up to 
j = k — 4, under the condition (a , e ) 7^ (0, 0). Now, we have a± — a 2 — ■ ■ ■ — a k -3 = 0. 
If we assume that a k -2 7^ then we can see from (|4.9| ) with j = k — 3 that ao = 0. Then 
the kth equation of (|4.10| ) reads 

8k(k + 2)a k _ 2 

k-l 60 = °' (4 - 20) 

which shows that e = since a k ^ 2 7^ 0. This contradicts the condition (b , a , d , e ) 7^ 
(0,0,0,0). Therefore, a k -2 must vanish. Next, let us assume that a k -i 7^ 0. Then we 
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can see from Q4.9| ) with j = k — 2 that ao = 0. Then the (k + l)-st equation of ( |4.10| ) 
becomes 

- 4(fc - l)a fc _ie = 0, (4.21) 

which shows that eo = since a k ~i 7^ 0. This contradicts the condition (6 , Oo> e o) 7^ 
(0,0,0,0). Therefore, must vanish. Then it has shown that a,\ = a% = ■ ■ ■ = 

a k -i = 0. Thus, we obtain 

V k = x k + a k y k . (4.22) 



(iii) When a 2 = (k(k — l)/2)ao and a = 1, the third and the fourth equations of (|4.10j ) 
yield 

M 32 g + M 31 b 

d ° = M33 ' (4 ' 23) 

1 U M M A3 M 32 \ f M 43 M 31 \ \ 

Then the condition (b , a , d , e ) 7^ (0,0,0,0) reads (&0) a o) 7^ (0,0). We obtain 
'k 



(sin^ ) J ± / ^k-j (cos^o) J 
(cosy? ) j ~ 2 (sinv9 ) J '" 2 



tan2^ = -— U) (4-25) 



from the recurrence relation (|4.9| ) under the condition (6 ,a-o) 7^ (0,0). Here, we take 
the positive sign for even k and the negative sign for odd k. The potential with ( 4.25f) 
is transformed into the separable form 

V k = (sec if ) k - V ± (cosec f ) k ~ 2 y k (4.26) 

by the rotation of coordinates defined by 

x — > xcos^o — l/sin^Oj ?/ — » x sin ^ + y cos y9 . (4.27) 

(iv) When a 2 = {?>k{2k — l)/2)a and a = 1, the second, the third, and the fourth 
equations of (|4.10| ) yield 



«o = -tt^o, (4-28) 

e » = { M « ^ Ml2 s| - Ml3 Ws ( Msl " Afe H) } V (4 ' 30) 

Then the condition (bo,ao,d ,eo) 7^ (0,0,0,0) reads 60 7^ 0. Let us put b = 1. Then, 
from the recurrence relation ( fOD , we can express a 4 , a 5 , . . . , ot k as polynomials of a 3 
whose coefficients are rational functions of k. Therefore, the rest of equations (|4.10|) , 
the number of which is 2k — 5, yield algebraic equations for a 3 . Now, the problem is 
whether the 2k — 5 algebraic equations for a 3 have common solutions or not. The two 
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algebraic equations of the lowest degrees, which are obtained from the fifth and the 
sixth equations of ( |4.10| ), are explicitly given by 

9(k - 2)k 2 {k + 2){2k - l)(5k - 2)G x {k) 



and 



+ 



4(fc - l)S(k) 2 
3{k + 2)(5fc - 2)(7fc - 6)G 2 (k)a 2 
4(fc - 2){k - l)(2k - l)(3fc - 2)R{k)S{k) 2 



2(fc + 2)(5fc-2)(7fc-6)G 3 (fc)at _ n f4 ^ 

(fc - 2) 3 (A; - l)k 2 (2k - l) 3 (3k - 2) 3 R(k)S(k) 2 1 ' ' 



9(k - 3)k(k + 2)(5fc - 2)G 4 (fc)a 3 



+ 



20(k - l)(3k - 2)R(k)S(k) 2 
3{k - 3){k + 2)(5k - 2){7k - 6)G 5 {k)a 3 3 
20{k - 2) 2 {k - l)k{2k - l) 2 {3k - 2) 2 R{k)S{kf 



(fc-3)(fc + 2)(5fc-2)(7fc-6)G 6 (fc)a| 
5(A;-2) 4 (A;-1)A; 3 (2A;-1) 4 (3A;-2) 4 J R(A;)5(A;) 2 ' 1 J 

where 

i2(jfe) = 74 - 151/t + 95fc 2 , S'(ifc) = -124 + 456A; - 519A; 2 + 241k 3 , 
G x (k) = 29952 - 260400A; + 932160A; 2 - 1727456A; 3 + 1706012A; 4 

- 810861A; 5 + 118529/fc 6 + 16870A; 7 , 

G 2 (k) = - 3170144 + 37193808A; - 190561760A; 2 + 553816792/c 3 

- 998755638A; 4 + 1151765545A; 5 - 843129587A; 6 
+ 37178381 lk 7 - 86588015/c 8 + 7446900A; 9 , 

G 3 (k) = - 1361952 + 20130192A; - 133944304A; 2 + 527295832A; 3 - 1360307178A; 4 
+ 2408313485A; 5 - 2978361002& 6 + 2565200757A; 7 - 1500811081A; 8 
+ 563283562A; 9 - 120201105A; 10 + 10736550A; 11 , 

Gi{k) = - 54723072 + 687806528& - 3817158400A; 2 + 12223132176A; 3 

- 24787677856A; 4 + 32872864764A; 5 - 28351948208A; 6 

+ 15180036291A; 7 - 4446138980A; 8 + 454376975A; 9 + 40055750A; 10 , 
G 5 (k) = 68988096 - 912699392A; + 5346987920A; 2 - 18143124368& 3 
+ 39260441572& 4 - 56288554120A; 5 + 53780522743& 6 

- 33457433369A; 7 + 12720106145/c 8 - 2562547775A; 9 
+ 186064500A; 10 , 

G 6 (k) = 56393856 - 913712256/c + 6725558432A; 2 - 29645084416/c 3 
+ 86929567784A; 4 - 178232080840 A; 5 + 261364852406A; 6 

- 275416174516A; 7 + 206176050353A; 8 - 106294583086& 9 
+ 35472298955A; 10 - 6761687100& 11 + 538285500& 12 . 
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The quantity called the resultant determines whether given two polynomials have a 
common root or not. Let us consider the two polynomials 

f(x) = a x n + a 1 x n - 1 + --- + a n (a ^ 0), (4.33) 
g(x) = b x m + b lX m - 1 + • • • + b m (b ^0). (4.34) 

The resultant of f(x) and g(x) is defined by the following determinant of order (m + n). 

Oq ai ■ ■ ■ a n 

ao ai ■ • • a n 



R(f,g) 



h ■■■ b m 

b h •■■ b m 



(4.35) 



b h b m 

Here, all the blanks represent zeros. It is known that the following theorem holds. 
Theorem 3 (see, for example, ||21||) The two polynomials f(x) and g(x) have a 



common root if and only if R(f,g) = 0. 

See |Appendix B| for the proof. The resultant of the two algebraic equations ( |4.31| ) and 
( |4.32| ) is computed to be fl22| 

289340 
~ 4096 

(k + 2f{k - 3) 4 (fc - 4) 2 (2fc - 3) 4 (3fc - 4) 4 (3fc - l) 4 (5fc - 6) 2 (5fc - 2) 29 (7A: - 6) 4 
X k 2 (k - l) 9 (k - 2) n (2A; - l) n (3A; - 2) 1G R(k) 6 S(k) 10 

x (29952 - 260400A; + 932160A; 2 - 1727456A; 3 

+ 1706012fc 4 - 810861A; 5 + 118529A; 6 + 16870fc 7 ) 
x (-2787216 + 23021920A; - 82567904A; 2 + 167598204A; 3 

- 211881739A; 4 + 173025983A; 5 - 91304549A; 6 

+ 30037885A; 7 - 5594600A; 8 + 450000fc 9 ) 2 , 



(4.36) 



which does not vanish for k > 5. Then it is concluded from theorem ^| that the two 
algebraic equations ( |4.31| ) and ([4.32|) do not have any common solutions. Therefore, 
there exist no common solutions among the 2k — 5 algebraic equations. This means that 
we have no solutions for a 2 = (3k(2k — l)/2)a («o — 1)- 

From the above arguments, we see that the case 5 considered here only yields the 
separable potential of the form 

V k = Ax k + By k , (4.37) 

i.e. there exist no new integrable potentials. 
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5. A list of integrable homogeneous polynomial potentials 



Although the computations in the previous section were performed under the implicit 
assumption that the degree of the potential, k, is greater than 4, we can carry out the 
computations for the cases k = 3 and k = 4 in the same manner. For k = 3, the 
expressions of My are given by 

(3a - a 2 )(45a - 2a 2 ) 



M 



ii 



M 21 = -30a 3 (3a - a 2 ), M 22 



112a 2 (3a;o — a 2 ) 



2 

M 4 i = -30a 2 a 3 , M 42 

8a 2 (3a;o — 16a 2 ) 



'33 



M :u = -9a a 2 - 7a 2 + 45a|, M 32 = 60a 2 a 3 , M 

-40a 2 , M 43 = 30a 3 (3a - 2a 2 



9(a - 2a 2 )(5a - a 2 ) 



Then the condition that the product MnM 22 M 33 M 4 4 vanishes gives the following 
relations between ojq and a 2 . 



45 



a 2 = 0, a 2 = 3a , a 2 = -^-«o, «2 = ^o, «2 
When a 2 = 0, we obtain the separable potential 



5a , 



«2 = T77«0- 

lb 



V 3 = a x 3 + a 3 |/ 3 



and when a 2 = 3a (a = 1), we obtain the potential 



V 3 = x 3 + 3xy 2 



«32T 



(5.1) 
(5.2) 
(5.3) 



which is transformed into the form of (|5.2|) by the rotation of coordinates defined by 
( |4.271 ). When a 2 = (45/2)ao («o = 1), we obtain an algebraic equation for a 3 , which 
is given by (|4.31| ) with k = 3, after the same computations as in the previous section. 
Then we obtain the potentials 



x 3 + 



45 
45 



x + -^%y + 



-xy 



nVTii 

2 ' 
27^^ • 



-y 



When a 2 



2 ° 2 
(l/2)a (a = 1), we obtain the potential 



x 



1 2 

2 Xy 



V3i 



-y 



When a 2 = 5ao («o = 1), we obtain the potential 



V 3 = x 3 + 5xy 



22^3 



(5.4) 
(5.5) 

(5.6) 

(5.7) 



The potentials (|5.5|), (|5.6j), ( p.7| ) are transformed into each other by proper rotations of 
coordinates. When a 2 = (3/16)«o («o = 1), we obtain the potential 



X 



— xy 
16 y 



(5.8) 
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The potentials ( |5.4| ), ( [5.8| ) are transformed into each other by proper rotations of 
coordinates. 

For k = 4, the expressions of M^ are given by 

(6a - a 2 )(42a - a 2 ) 



M 



11 



M 2 i = -lla 3 (6a - oc 2 ), M 22 
M 3 i 



3 

80a 2 (6a;o — a 2 ) 



3 

540a a 2 — 16af + 99^ — 1512a «4 + 228a 2 a: 4 



6 

M 32 = -80a 3 (3a - a 2 ), M 33 = 80a a2, 

* 8(32a§ - 27a 2 , - 24a 2 a 4 ) 
M 41 = -2a 3 (6a + 37a 2 - 48a 4 ), M 42 = ^ 2 - -2 

O 

n 64a 2 (3a — 4a 2 ) 
M 43 = 24a 3 (7a -a 2 ), M 44 = 2 ^ 



3 

Then the condition that the product MnM 22 M 33 M 44 vanishes gives the following 
relations between a and a 2 . 

3 

«o = 0, a 2 = 0, a 2 = 6«0) a 2 — 42a , a 2 = ^ a o- (5-9) 
When a , we obtain the potential 

V 4 = y\ (5.10) 
and when a 2 = (a^ = 1), we obtain the potential 

V 4 = x 4 + a 4 y A . (5.11) 
When a 2 = 6a (ao = 1), we obtain the potentials 

V 4 = x A + 6xV + a 3 xy 3 + % 4 , (5.12) 

16 

V 4 = x 4 + 6xV + 8y\ (5.13) 

The potential ( 5.12j ) is transformed into the form of ( p. 11 ) by the rotation of coordinates 
( [4.271 ). When a 2 = 42a («o = 1) ; we obtain three algebraic equations for a 3 , two of 
which are given by ( [4.3 1| ) and ( |4.32| ) with k = 4, after the same computations as in the 
previous section. They have common solutions, which is indicated by the fact that the 
resultant ( f4.36| ) vanishes when k = 4. Then we obtain the potential 

V 4 = x 4 + 42xV + 28Vl0ixy 3 - 48y 4 . (5.14) 
When a 2 = (3/4)ao («o = 1)) we obtain the potential 

V, = x 4 + \x 2 y 2 + l -y\ (5.15) 

The potentials (|5.13|) , (|5.14| ), Q5.15| ) are transformed into one another by proper rotations 
of coordinates. 

As a consequence, we obtain the following list of integrable 2D homogeneous 
polynomial potentials with a polynomial first integral of order at most 4 in the momenta. 
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with a polynomial first integral linear in the momenta 

V k = (x 2 + y 2 ) k/2 , k = even, 
with a polynomial first integral quadratic in the momenta 



Vk 



r + x 



r — x 



V k = Ax" + By* 



with a polynomial first integral quartic in the momenta 



lb 



V* 



y/3't 



x + -xy + 



-r 



V4 



4 1 ^ 2 2 1 1 4 
x + -a: y + -y . 



(5.16) 



(5.17) 



(5.18) 



As far as the authors know, no one has discovered any polynomial first integral which is 
genuinely quintic or higher orders in the momenta. It is still an open problem whether 
or not there exist such polynomial first integrals. 
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Appendix A. Proof of the properties 1 and 2 

Appendix A.l. Time reflection symmetry - proof of the property 1 
The system ( |1 . 1| ) is invariant by the time reflection 

t^-t, x^x, y^y, p x -> -p x , py^-py. (A.l) 

If $(x,y,p x ,p y ) is a first integral of the system (Oj), then $(x, y, — p x , — p y ) is also a 



first integral because of the time reflection symmetry of the system. Note here that 
every first integral $ can be decomposed into the sum of $ ev en and $ dd, given by 

<P(x,y } p X} p y ) + $(x,y,-p x ,-py) , 

^even — ~ 

and 

$(x, y, p x , Py) - <2>(x, y, -p x , -p y ) 

*odd — • l A -<jJ 

We can see that ^ ev en is a first integral which is even in the momenta and that $ dd 
is a first integral which is odd in the momenta. That is, if $ = $ G ven + ^odd is a first 
integral of the system ([Lip, then $ e ven and $ Q dd are also first integrals. Therefore, we 
can assume that a first integral is either even or odd in the momenta from the beginning. 
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The system ( |1 . 1| ) with a homogeneous polynomial potential of degree k is invariant by 
the scale transformation 



t 



<T~H, 



X 



a 2 '^x, 



y^a^'%, Px 
In general, a system of differential equations 

^ = F i (x 1 ,x 2 , ■ ■ ■ ,x n ), (i = l,2 

is called a scale invariant system if it is invariant by the scale transformation 

t-KT~H, Xi^a 9l Xi, (i = 1,2,..., N) 

with an arbitrary parameter a and proper constants A function $(xi,x 2 , ■ 
said to be weighted-homogeneous with a weight M if it satisfies 

$(a 9l x u a 92 x 2 , . . . , a 9N x N ) = cr M ${x 1 ,x 2 , . . . , x N ). 

Suppose now that the scale invariant system ( |A.5|) has a polynomial first integral 
which can be written in the form 

in 

where each polynomial <I? m is weighted-homogeneous with a weight m. The scale 
transformation (|A.6|) transforms the first integral (|A.8| ) into 



o- k/ik - 2) Px , 



N) 



Py -> <T k/ik -% y . (A.4) 

(A.5) 



(A.6) 

,x N ) is 

(A.7) 



m 



(A.9) 



which is again a first integral for an arbitrary a because of the scale invariance of the 
system. Then, it is concluded that each polynomial $ m is a first integral. Therefore, we 
can assume that a polynomial first integral is weighted-homogeneous from the beginning. 



Appendix B. Resultant — proof of theorem || 

The resultant is an algebraic tool for elimination of a variable between two algebraic 
equations and gives the condition that the two algebraic equations have a common root 
(theorem |3|). See for more details of the resultant and its applications. 



Proof of theorem ^ Suppose that the two polynomials have a common root, say a. 
Then the following simultaneous algebraic equations hold. 



m ~ 2 /(«) 



a 



a a m+n - 2 + ai a m+n - 3 + 



+ a n a 
+ a n a 



m—l 



m-2 








/(«) 



a a n + a\a n 1 + 



a" 1 g(a) 
a n ~ 2 g(a) 



9{a) 



+ a n = 

b a m+n - 1 + M m+ "~ 2 + • • • + b m a n ~ l 



(B.l) 



boa 



m+n—2 



+ b\OL 



m+n—3 



H h b m a 



n-2 



b a m + M m_1 + ••■ + & 
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If we multiply the Ith column of the resultant by a 1 and add them to the (m + n)-th 
column, then all the elements of the (m + n)-th column of the resultant vanish because 
of (Pip . Therefore, R(f,g) = 0. 

Let us next assume that R(f, g) = 0. Then the m+n row vectors of the resultant are 
linear dependent. Let at, a 2) . . . , a m , bx, b 2 , . . . , b n be the row vectors of the resultant. 
Then the relation 

m n 

J2cia i + J2 d J b 3 = ( B - 2 ) 

i=l j=l 

holds with (ci, . . . , d n ) ^ (0, . . . , 0). Multiplying the Zth element of ( |B.2|) by x m+n ~ l and 
adding them up, we have 

m n 

c lX m - l f(x) + djx n -ig(x) = 0, (B.3) 

i=l j=l 

which is reduced to 

h(x)f(x) = -k(x)g(x), (B.4) 

where Y^iL i CiX m ~ l = h(x),Y^j=idjX n ~ :i = k(x). Let deg/(x) denote the degree of a 
polynomial f(x). Then, we have the relation 

deg k(x) < n — 1 < n = deg f(x). (B-5) 

If f(x) and g(x) have no common factors, then f(x) must be a factor of k(x) because 
of ( |B.4| ). This contradicts equation ( |B.5| ). Therefore, f(x) and g(x) have at least one 
common factor, i.e. they have at least one common root. This completes the proof of 
theorem |^. 
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